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^ ■ Abstract 
> 

^ ! We consider an astrophysical system with a population of sources and a 



population of lenses. For each pair of source and lens, there is a thin on-axis 
tube-like volume behind the lens in which the radiation flux from the source is 
greatly increased due to gravitational lensing. Any objects (such as dust grains) 
which pass through such a thin tube will experience strong bursts of radiation, 
i.e.. Extreme Gravitational Lensing Events (EGLE). We study the physics and 
statistics of EGLE for the case in which the shear is larger or comparable to the 
finite source size. We find that the presence of shear has only a small effect on 
the EGLE statistics. 
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1. Introduction 

Gravitational lensing in the context of direct observations of lensed sources has been 
studied extensively in past years (piandford and Narayan 1992| , Schneider et al. 1992 



and references therein). Recently, we proposed a new way of looking at gravitational 
lensing, the statistics of extreme gravitational lensing events (EGLE), which may have 
significant astrophysical effects ( [Wang and Turner 1996| , hereafter referred to as Paper 1). 



We introduced the basics of EGLE in Paper 1, where we considered the case in which the 
finite source size is more important than shear. In this paper, we discuss the case in which 
shear can not be neglected. 

Extreme gravitational lensing events (EGLE) occur when the source, lens, and observer 
(or target) are nearly on-axis. In the statistics of EGLE, we consider an astrophysical 
system of diameter Dc, in which both the sources and lenses are uniformly distributed, and 
the targets can be the fragile objects in the same system, such as dust grains, etc. 

If the target moves a distance d away from the line connecting the source and the lens, 
it is equivalent to the source moving an angular distance of y from the optical axis (the line 
connecting the lens and the target). Measuring y in units of the angular Einstein radius, 
we have 

where Djs, -Dg, and are angular diameter distances between the lens and source, target 



and source, target and lens respectively. 6'e = \j2Ra,D^f,/ {D^D^) is the angular Einstein 
radius. -Rs = 2GM is the Schwarzschild radius of the lens with mass M. Eq.(0) is the small 
angle approximation, valid for d ^ D^. The dimensionless radius of a source with physical 
radius p is defined to be = p/ {Ds9e)- 

For a given pair of lens and source, the shear 7 due to other mass inhomogeneities near 
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the line-of-sight is the same order of magnitude as the optical depth for microlensing, r, the 
probabihty that the source is lensed. We find [see Appendix] 



C (x) is the Riemann zeta function. Not surprisingly, the statistics of EGLE is much more 
complicated for sheared lenses than for isolated lenses (7 = 0). 

The lens model for a perturbed Schwarzschild lens (i.e., point mass lens with shear) 
has been investigated by Chang and Refsdal (1979, 1984) and Subramanian and Chitre 
(1985). The magnification of extended sources by a perturbed Schwarzschild lens has been 
studied by Chang (1984), Schneider and Weiss (1987). In Paper I, we showed that shear is 
not important for 7 < i? <^ 1. 

In this paper, we study the statistics of EGLE for sources with small dimensionless 
radius R and slightly sheared lenses (7 <^ 1), with 7 > i?. We generally follow the notation 
and conventions of Schneider et al. (1992) and Paper I. 



In this section, we study the basic properties of a perturbed Schwarzschild lens. In 
the absence of shear, the caustic (which corresponds to infinite magnification) of a point 
mass lens in the source plane is a point (intersection with the optical axis) . If the shear of 
a macrolens at the location of the point mass lens is < 7 < 1, the caustic in the source 
plane changes from a point to an astroid shaped closed curve with four cusps, assuming 
that the surface mass density of the macrolens is zero at the location of the point mass lens. 

Let us define 




(2) 



2. The perturbed Schwarzschild lens 



Y = 



y 



1-7 



(3) 



1 + 7 
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The caustic curve intersects the Yi axis at ±(1 — A), the Y2 axis at ±(1 — A)/-\/A. For 
7 <S 1, the caustic is described by 

Yi = -27 cos^ 0, Y2 = 27 sin^ 0. (4) 

As we discussed before, moving the observer away from the hne connecting the source and 
lens is equivalent to moving the source relative to the optical axis. When a source moves 
off the optical axis, if it moves along Yi = or F2 = 0, it crosses a cusp; if it moves along 
Yi = Y2, it crosses a fold. All observables for a source moving off the optical axis in an 
arbitrary direction fall between the two limits set by Y1Y2 — and Yi — Y2. We therefore 
restrict our attention to these two relatively simple cases below. 

The lens mapping equation can be solved exactly for Y1Y2 — 0. (Chang & Refsdal 
1979 & 1984, Subramanian & Chitre 1985) For a point source moving along the Yi axis, its 
magnification is 

1 {1 - A){1 + A) - Y,^ 



/Xp(|Fi|<l-A) 



;i + 7)2 A[(i-A)2-y;: 



21 ' 



(1 + 7)' a^i^^Ima [r^ _ (1 _ A)2] 

For a point source moving along the Y2 axis, its magnification is 

^ 1-A\ 1 (1-A)(l + A)+Ar22 



VA ; (1 + 7)2 A[{l-Ay-AYi] ' 



„ 1-A\ 1 Y2[Yi + {3-A)] 

^^V Va) + ^W^4[AYi-il-Ar] ^ ^ 

For 7 <^ 1, 1 — A ~ 27, Y ~ y. Let us define 

^ 27' 27 ^ ' 



Then, for ^2 = 0, 



, ^ 1 1 , ^ 1 1 + 27%? , , 



71-Fr -p^'-'-^^ 72(^2-1)^17- 



2y? 
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For a small source with R < 0.2, the maximum magnification is given by j2p {\y^\ < 1), with 

- \ / — \0.73 

(1 — Hi) ( 1.2i?j . For yi = 0, just replace i/i with ?/2- Figure 1 shows the lightcurves 
for R = 0.01, 0.1, 0.4, 1 and y2 = 0. Note that the smaller the source size compared to 
the shear, the larger the magnification of the source near the caustic; the peaks in the 
lightcurves correspond to caustic-crossing. 

For a source moving along the diagonal line yi = y2, (7 — y) gives the source's distance 
to the caustic. For a small source with R < 0.2, the maximum magnification is given by 
the point source magnification /ip(y), with (0.5 — y)= 0.67R. The magnification of a point 
source near the caustic is given by ( |(Jhang and Kefsdal 198'^ ) 

/i - 0.5) ^ i -=^. (9) 
7 VO.5 - y 



For a very small source (i? < 0.02) ([Chang and Refsdal 1984|) , 
X 1 0.3 , y-0.5 

if \xo\ < 1, Zmax = -Xq] elsC Z^ax = 1- 



2 

TT 



Z^ 



-1 \^-{z + xo)' 

Figure 2 shows the lightcurves for R = 0.01, 0.1, 0.4, 1 and y2 = yi- Again, the peaks in 
the lightcurves correspond to caustic-crossing. 

For R > 0.01, we fit /Umax {Ri l) with the following formulas: 



7/^11 



l-w{R) w{R) . , 

^ + —^m-, cusp-crossmg (2/2 = 0); 

R 2R' 



w{R) 



R Vl 



0.3555 

w{R), fold-crossing {y2 = Vi)- (H) 



w = exp 



—R{1 + R)/2 is the weight function. The above fit formulas underestimate the 
maximum magnification by about 20-30% when the size of the source is comparable to the 
size of the region enclosed by the caustic curve. Note that for 7 = (i? = 00), we recover 
A^max = 2/-R, the exact result used in Paper I. For R < 1 [i.e., R < 27], /imax < 2/R. 
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3. EGLE volume statistics: general description 



Let us consider a mixed population of uniformly distributed sources (with number 
density ns) and lenses (with number density wl) confined in a volume of diameter D^. We 
are interested in computing the volume fraction of space in which the flux exceeds / due to 
gravitational lensing, i.e., the volume fraction of space in EGLE tubes with flux > /. 



R — Rn 



t 



x{t + x) ' 



7 



Note that in EGLE statistics, both 7 and R depend on distances. We write 

jotx, t + X < 1, 
'yotx/{t + x), t + x>l. 
where 

P 



(12) 



Rn 



7o = 3.7to. 



p is the physical radius of the source, Tq is the optical depth at Ds — D^. 
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For a source with luminosity Ls, the unlensed flux /o — Ls/I^ttD"^]. As in Paper I, we 
deflne 



a 



PmJxi-R) 


2 




\ iVRo) 1 


[ f/fo \ 






. (///min) _ 



(14) 



p2 {AttD^JJ ' 

where /min = -^s/[47r£)^]. a measures the maximum magniflcation of the source relative to 
the flux /, in the limit of zero shear. 

We use the simple but reasonable approximation: 

max 



i^e{y,R,i) - ' 



(15) 



Afmax elsewhere 

where //p(y,7) is the point source magniflcation in the presence of small shear, and //max is 
the maximum magniflcation. For 7 = 0, we recover the approximation used in Paper I. Let 



7 



us write [see Eq.(|TT])] 



/^max 



(16) 



w (-R/7) is a weight function. For 7 <^ i?, we recover the exact formula used in Paper I. 

For a given pair of lens and source, we define the EGLE tube with minimum flux / 
to be the small tube-shaped volume behind the lens in which the flux exceeds / due to 
gravitational lensing. For a given pair of lens and source with separation x = D^s/Dc, 
the length of the EGLE tube is tmax = D^d^'^/Dc, i.e., the EGLE tube extends from t = 
to t = tmax- The minimum magnification needed at a point inside the EGLE tube is 



= f/fo ={t + xy ///min, the unlensed flux /o = Ls/[AttD^] = /mi„/(t + xf. The length 



of the EGLE tube is given by setting n = /imax(^, x). 

The cross-section of the EGLE tube at a distance of behind the lens along the line 
SL which connects the source and the lens, a{Dd,f), is the area perpendicular to SL in 
which the magnification exceeds = ///o, i.e., the flux exceeds /. The cross-section of the 
EGLE tube can be written in the form 



a = a/[d{fi)/y{fi)]'^ is the cross-section of the EGLE tube in units of Einstein radius squared 
[see Eq.(|l])]. For = 7 = 0, a = vr(i^(/i), jj, ~ l/y, hence a* = 1. For 7 7^ 0, a* = o-*{-fij,). 

The volume fraction of space occupied by EGLE tubes in which the flux exceeds /, JFl, 
can be computed using the same method as in Paper I. 

For a given lens and source pair, the associated EGLE tube with minimum flux / has 
the volume 





(17) 




2 



1 ft' 
X Jo 



t 



dDd(r{Dd,f) = 2nRsDl 



dt 



it + x)3 



a*{t,x). 



(18) 
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The volume fraction of space in EGLE tubes of flux > / for a given lens is 

Vdf) = 47rns ^ d^ds Dl Vsl(/). (19) 
Jo 

The volume fraction of space in EGLE tubes of flux > / is 

Mf) = nMf) = J^df,R = l = 0) I{a, i?o/7o), (20) 

where J-'hif, -R = 7 = 0) is the volume fraction of space in EGLE tubes of flux > / for point 
source in the absence of shear (see Paper I), 

^L(/,i? = 7 = 0) = 47r2nsnLi?sI^,^(^^^ , (21) 

and 

/(a,/?o/7o) = ^.^'^^'V^ , = 2 f'dxx f^^^dt-^^a*it,x). (22) 
^l(/,P = = 7) Jo Jo [t + xy 

I{a,Ro/'j()) is the modification factor in EGLE volume fraction due to non-zero shear 
or/and finite source size. For 70 = 0, I{a, Rq/'Jo) = I{o:) (see Paper I). 

For given finite ratio -Ro/70, -^(«5-Ro/7o) has the following asymptotic behavior. At 
very large a, both R and 7 become negligible, and J(a, i?o/7o) approaches 1. If a is not too 
small, EGLE tubes with flux > / exist for arbitrary source- lens separations (0 < D^s < -De). 
When a is sufficiently small, EGLE tubes with flux > / do not exist for source-lens 
separations larger than the maximum D"^!^^ < D^. For decreasing a, the EGLE tubes with 
flux > / decrease in length and they require smaller source-lens separations; when the 
EGLE tubes are sufficiently short and require sufficiently small source-lens separations, 
i? > 7, the finite size of the source dominates [see Eq.([T^)], and /(a, -Ro/70) approaches 
/(a). 
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4. EGLE volume statistics for 7 7^ 0, i? 7^ 

In the presence of shear, the EGLE tube has very comphcated topology. For given R, 
7 7^ always leads to smaller /imax [see Eq.([TT|)], hence shorter EGLE tubes. The cross 
section of an EGLE tube is determined by the shape of the caustic curve [with both cusps 
and folds] in the source plane. For 7 7^ 0, the volume in which the flux produced by EGLE 
exceeds the critical value has the topology of a filled tube over most of its length, but near 
the end of the tube a complex geometry develops refiecting the cusp and fold natures of the 
caustic. 

For a point source, the EGLE tube has an astroid shaped cross section until = D^, 
with given by ///o(/^d) = 1/7. For > D'^, the EGLE tube is hollow at the center, 
with walls tapering off to infinity away from the center; the walls are thicker and rounded 
in the cusp directions, and thinner with steep outside surfaces half-way between the cusp 
directions. The walls taper off more quickly half-way between the cusp directions than in 
the cusp directions. 

For a finite source, if the source is sufficiently large {R > 27), the EGLE tube ends 
without the shear-induced hollow and the effect of shear is negligible. If i? < 27, the EGLE 
tube has a topology similar to the R = case, except that the shear-induced hollow has 
walls which end at a finite distance behind the lens, the walls stretch furthest in the cusp 
directions, and shortest half-way between the cusps; i.e., the end of the EGLE tube is a 
hollow volume, with thinning walls which terminate in four rounded legs centered in the 
cusp directions. 

The volume of an EGLE tube for 1^7>i?>0is between the volume of an EGLE 
tube assuming cusp-crossing only and the volume of an EGLE tube assuming fold-crossing 
only. Compared to the point source magnification in the absence of shear (dashed curve in 
Figs. 1-2), for a given magnification of the source, the sheared lens (with 7 > i?) leads to 



- 10 - 



a much larger cross-section for cusp-crossing [see Figure 1] and a smaller cross-section for 
fold-crossing [see Figure 2]. Therefore, the volume modification factor /(a,_Ro/7o) should 
be smaller in the fold-crossing case than in the zero shear case; but it can be larger in the 
cusp-crossing case than in the zero shear case, because of the significant thickening of the 
EGLE tubes at the ends away from the lenses. 

In this section, we use Eq.(0) to compute the EGLE tube lengths, and 
f^piy) = fJ'pi.V: 7 7^ 0, i? = 0) to compute the EGLE tube cross-sections. 

The cross-section of the EGLE tube can be expressed in terms of the dimensionless 
parameter a* = fi'^a/n [see Eq.(|l^], where a = a/[d{fj.)/y{fj,)]'^ is the cross-section of the 
EGLE tube in units of Einstein radius squared. 



1/(2) and are roots of fie{y,'y) = ^, y(2) > > (if fie{y,l) = only has one root, then 
we denote it by ?/(2), and set = 0). For 7 <^ 1, 7/ie(?/,7) only depends ony = y/{2j). 
Then 



For < 0.1/7, the shear on the lens has negligible effect on the lightcurve, hence 
(T* = 1 as in the i? = 7 = case. For fi > O.I/7, we can compute a* using the point 
source and small shear lightcurves shown in Figs. 1-2, for cusp-crossing and fold-crossing 
respectively. For cusp-crossing, we can compute a* analytically using Eq.(H). We find 



where 6 = 7/i. Note that for 5 = (i.e., 7 = 0), we recover a* = 1. For fold-crossing, we 
have to find the roots of 7 fie{y, 7) = Se{y) = ^ numerically to compute a*. 




(23) 





6 < 1 



(25) 
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To estimate the effect of shear on the EGLE volume fractions, we write 



where A = 57r/32 is 



Iavg{Ro, 7o) 




(26) 



(270)^]. The averaged EGLE volume fraction is only slightly larger than the EGLE volume 
fraction assuming fold-crossing only, because most of the caustic is fold. 

In Fig. 3, we show the modification factor in EGLE volume fraction due to finite 
source size and non-zero shear, Iavg{a, Ro/jo), for -Ro/70 = 10, 1, 0.1, 0.01, 0.001. The 
R0/10 = 10, 1 curves are indistinguishable from the 7 = curve (dashed line). Note that all 
curves in Fig. 3 converge to the 7 = curve at small and large a, as discussed at the end of 
last section. Note that the effect of 7 7^ decreases at both extremes of large a [small flux 
or high magnification] and small a [high flux or low magnification]. Clearly, the presence of 
shear has only a moderate effect on the EGLE volume fractions. 

For objects which enter an EGLE tube perpendicular to the lens-source line, the EGLE 
durations are longest in the cusp directions, and shortest half-way between the cusps. Since 
the EGLE durations also depend on the angle at which an object enters the EGLE tube, 
the effect of shear on the EGLE durations should be negligible statistically also. 



We find that the presence of small shear [comparable or larger than the dimensionless 
source size] has only a small effect on the EGLE volume fractions. Specifically, the EGLE 
volumes are typically decreased by factors of a few for 70 ~ 10^-Ro [the characteristic 



[which measures the maximum magnification of the source in the limit of zero shear relative 
to the flux /] and otherwise unaffected, see Fig. 3. This means that the 7 = results given 



5. Summary 



dimensionless shear 70 and source size Rq are defined in Eq.(P^] at moderate values of a 
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in Paper I will be adequate for most order of magnitude considerations. 

It is quite reasonable that the effect of < 7 <^ 1 is modest since the convergence 
provided by the EGLE lens itself constrains the total flux concentrated into the EGLE 
tube. The small shear considered in this paper only rearranges this flux somewhat into a 
slightly modified EGLE volume, especially near the end of the EGLE tube. 

Y.W. is supported by the DOE and NASA under Grant NAG5-2788. E.L.T. gratefully 
acknowledges support from NSF grant AST94- 19400. 



Appendix: Derivation of 7(r) 

Here we derive the shear 7 due to other mass inhomogeneities near the line-of-sight 
for a given pair of lens and source, in terms of the optical depth for microlensing, r, the 
probability that the source is lensed. 

We assume that the EGLE lens [the lens under consideration] is embedded in a 
distribution of identical lenses with surface number density S, each lens has the same 
angular Einstein radius 9e associated with it. The total shear on the EGLE lens is 



where A = l/y/nT, is the dimensionless distance between two adjacent lenses. Using 

T = 7rE6'|, and E^i J'^^^ = C(3/2) ~ 2.6124 [where C{x) is the Riemann zeta function], we 




(1) 




obtain 




(2) 
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Fig. 1. — The lightcurves for y2 = (cusp-crossing), R/{2^) = (dotted line), 0.01, 0.1, 0.4, 
1 (solid lines). The dashed hne is the lightcurve for i? = 7 = 0. 

Fig. 2.— The lightcurves for = Vi (fold-crossing), R/{2-f) = (dotted line), 0.01, 0.1, 0.4, 
1 (sohd lines). The dashed hne is the lightcurve for = 7 = 0. 

Fig. 3. — The modification factor in EGLE volume fraction due to finite source size and 
non-zero shear, Iavg{ot-i -Ro/7o)/-^(«, Ro/jo — 00), for i?o/7o — 10 (solid line), 1 (dotted line), 
0.1 (short dashed line), 0.01 (long dashed line), 0.001 (dot-short dashed line). 
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Fig. 1. — The lightcurves for y2 = (cusp-crossing), R/{2^) = (dotted line), 0.01, 0.1, 0.4, 
1 (solid lines). The dashed hne is the lightcurve for = 7 = 0. 
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Fig. 2.— The lightcurves for = Vi (fold-crossing), R/{2-f) = (dotted line), 0.01, 0.1, 0.4, 
1 (solid lines). The dashed hne is the lightcurve for = 7 = 0. 
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Fig. 3. — The modification factor in EGLE volume fraction due to finite source size and 
non-zero shear, Iavg{a, Ro/^o)/ I{a., -Ro/to — oo), for Ro/^o — 10 (sohd hne), 1 (dotted hne), 
0.1 (short dashed line), 0.01 (long dashed line), 0.001 (dot-short dashed line). 



